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ABSTRACT 

We analyze the weak lensing data of the VIRMOS imaging survey using pro- 
jections (called E and S-modes) of the two independents observed correlation 
functions. The E-mode contains all the lensing signal, while noise and system- 
atics contribute equally to the E and B modes provided that intrinsic alignment 
is negligible. The mode separation allows a measurement of the signal with a 
\/2 smaller error bars, and a separate channel to test for systematic errors. We 
apply various transformations, including a spherical harmonic space power spec- 
trum Cf and Cf , which provides a direct measurement of the projected dark 
matter distribution for 500 < / < 10 4 . 

Subject headings: cosmology: dark matter - lensing 



1 Based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which is operated by 
the National Research Council of Canada (NRCC), the Institut des Sciences de l'Univers (INSU) of the 
Centre National de la Recherche Scientifique (CNRS) and the University of Hawaii (UH) 
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1. 



Introduction 



The measurements of the alignment of distant galaxies from ground based data (Van 
Waerbeke et al. (2000); Wittman et al. (2000); Bacon et al. (2000); Kaiser et al. (2000); Maoli 
et al. (2001)) and from space (Rhodes et al. (2001)) provided the first convincing evidence 
of gravitational lensing by large scale structures. In a more recent work, Van Waerbeke 
et al. (2001) measured the cosmic shear signal using different statistics, and shown the 
remarkable agreement between them, which demonstrated the gravitational lensing origin 
of the signal. However an indication of a remaining systematic were found in the aperture 
mass measurements. It could be either due to an imperfect Point Spread Function correction, 
and/or to an intrinsic alignment of galaxies arising from spin-spin correlation for instance 
(Crittenden et al. (2000)). Although a robust and ultimate method to clean the measured 
lensing signal from contamination would be the measurement of the cross-correlation between 
different source planes, it is in principle possible to clean existing data using a curl- free/curl 
modes decomposition (called E and B modes). The reason is that gravitational lensing 
produces curl-free shear patterns only, while systematics and intrinsic alignments contribute 
to the B mode as well (Crittenden et al. (2000)). This paper is a tentative to measure these 
modes separately in the lensing signal measured in the VIRMOS-Descart 2 weak lensing 
survey (Van Waerbeke et al. (2001)) and to establish the lensing origin of the signal on a 
quantitative base. 

The first Section establishes some basic relations, and defines the mode decomposition 
of the shear correlation function. The second Section shows the mode measurement on inte- 
grated shear correlation functions which are trivially related to windowed variances measured 
by previous authors. The last Section shows the power spectrum measurement for the two 
modes. The theoretical background of this work was essentially developed in Crittenden 
et al. (2000) and will not be detailed here. 



Weak lensing shear is a spin 2 polarization field, and can be described by a traceless 
two by two matrix at every point, 
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Correlation Functions 




(1) 



2 http://terapix. iap.fr/Descart/ 
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This matrix has two degrees of freedom at every point, much like a vector field. In analogy 
to a vector field, there are two coordinate invariant projections, a divergence and a curl. The 
former is called E-mode in analogy to the electric field, and the latter called I?-mode which 
is divergence free and analogous to a magnetic field (Zaldarriaga & Seljak (1998); Seljak & 
Zaldarriaga (1998)) 

The weak lensing shear can be expressed in terms of the mean expectation of source 
ellipticities and alignments 

7l = (ecos(20)), 7 2 = (esin(2#)> (2) 

where 9 is the angle between the major axis of the source galaxy and the x axis, and 
e = (a — b)/(a + b) is determined by the major axis length a and minor axis length b. 
Following Miralda-Escude (1991); Kaiser (1992) we define the shear components ( 7 t, 7r ) in 
the frame of the line connecting a pair of galaxies. In equation (2) this corresponds to 
defining the galaxy alignment 6 in the frame of the pair separation. 

To date, all shear correlation function analysis have been performed from measurements 
of ( r yt(x) r Yt(x + r )) and ( 7r (x) 7r (a; + r)). The spin-2 correlation function can then be decom- 
posed into two coordinate invariant components, £ + and 

f+(r) = ^(lij(xhij(x + r)> = ( 7t (x) 7 i(x + r)> + (-f r (x)j r (x + r)). (3) 

£-0*0 = {lt{x)lt(x + r)) - ( 7r (x) 7r (x + r)>. (4) 

To test for systematic errors in weak lensing mass reconstruction, one usually rotates all 
galaxies by 45 degrees, upon which the weak lensing signal should disappear. Under this 
rotation (3) is unchanged, while (4) changes sign. 

Weak lensing arises from a gravitational potential. The statistics of this scalar field can 
be described by a single correlation function, which means there must exist a degeneracy 
between the two correlators (3,4). Crittenden et al. (2000) have shown how to transform 
these two correlators into a pure .E-mode which contains all the lensing signal, and a pure 
-B-mode which should contain only noise. Any remaining systematic effects is expected to 
contribute to both E and 5-mode: residuals in the Point Spread Function corrections will 
contribute equally, while the intrinsic alignment correlation of the galaxies is expected to 
give a higher amplitude in E than in B. How much higher is still debated in the literature 
(Crittenden et al. (2000),Croft k Metzler (2000), Catelan et al. (2001), Heavens et al. (2000)). 
However, a general agreement among the literature is that our survey is deep enough so that 
intrinsic alignment has a negligible contribution (less than 10%, see for instance Pen et al. 
(2000)). 
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Following Crittenden et al. (2000), we define 



£'(r)=£_(r)+4 




(5) 



The E and -B-type correlators are given as 
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The sum is local, ^ E + £ B = £+, while the difference depends on the difference correlator 
at larger radii. Rotating all galaxies by 45 degrees swaps £ E and £ B , but leaves the sum 
unchanged. If we rotate only one galaxy in each pair by 45 degrees, both E and B correlators 
should disappear. In practice, we only know the correlators out to a finite radius r\. The 
integral in (5) which should go to infinity must thus be truncated at r\. The integral from 
r\ to infinity results in two integration constants which affect the correlation function at 
smaller radii. The integration constant for the last term in (5) is weighted by (r/ri) 2 , and 
should not affect short range correlations significantly. The middle term results in a straight 
constant £, which can be added to the E correlator and subtracted from the B correlator. 
The decomposition (6) is degenerate under a change £ E (r) — > £ E (r) +£, £ s (r) — > £ B (r) — £. 
This degeneracy traces back to the underlying map decomposition degeneracy. A constant 
shear signal generates a constant correlation function. But this constant shear is zero under 
differentiation, and cannot be classified unambiguously as curl or div type. In practice, we 
integrate as far as we have data, and parameterize all lack of knowledge in terms of the 
integration constant £. If one wishes to test the null hypothesis that there are no B modes 
as predicted by weak lensing, one can marginalize the results of all possible integration 
constants. If B cannot be set to zero at all scales for any value of the integration constant, we 
reject the hypothesis that there are no B modes. This would be an indicator of instrumental 
systematics. Since there are two integration constants, one should marginalize over both. 
One of the constants only affects large scales, so at smaller scales only £ is of importance. 

The raw correlation functions £ + , £_ are measured by averaging all pairs of galaxies at 
a given separation r. In practice, we binned the separations into intervals of three pixels, 
corresponding to about 0.6 arc seconds, for a total of 10000 bins corresponding to 1 2/3 
degrees separation. Each pair is given a statistical weight defined by: 



where Oi and <jj are the ellipticity measurement r.m.s. of the galaxies i and j. The quantity 
a e is the r.m.s ellipticity estimated over all the galaxies. This is different from the weighting 
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scheme expressed in Eq.(7) in Van Waerbeke et al. (2000) but it gives similar results. The 
sum of weights of pairs in each bin provides a total weight w for each bin. 

A crucial step is now the calculation of the noise covariance matrix. Since the number 
of pairs is the square of the number of actual galaxies, one might expect the statistical error 
of each correlation bin to be slightly correlated with other bins. In practice, we have not 
been able to find such a correlation (because the bins are small and the noise dominates 
over the signal), we thus treat each correlation bin as uncorrelated with the others. Each 
of the correlators now has mutually identical, independent and radially uncorrelated 

bins. The mapping (5) correlates bins at different separation, while the transformation (6) 
now correlates the E and B correlators. We compute the noise covariance matrix as follows. 
The variances are discretized as a one dimensional array e = {a 2 (r,j)}. We will use the 
subscript i to index this one dimensional bin array. Since we used 10000 bins Vi to compute 
the correlation functions, we can define the 20000 elements vector & = (£ + (rj), £_(rj)). The 
noise covariance matrix = — — of the raw correlators is a 20000 2 matrix 

which is e on the diagonal and zero elsewhere. The transformation into E and B correlators 
from (6) is a linear operation on the local correlators £j. We again define the 20000 elements 
vector £f B = (£ B (rj), £ B (rj)). Equation (6) implicitly relates them through a transformation 

CP = T,r r (8) 

Since each bin in £ EB is extremely noisy, we have re-binned it into 7 logarithmic intervals 
each a factor of two wide, denoted This re-binning can be represented as a 14 by 20000 
projection operator £j£ = P k i£ BB . The binned noise covariance matrix N b is 

Nl = <(£? - - (&»> = PuT l3 N jk T{ Pt m - (9) 

We show the result of the two projections in figure 1. The error bars are the square root 
of the diagonal entries of N b (9), and correspond to one a error bars keeping all other data 
points fixed. We have normalized the degeneracy £ such that ^ B is zero at the bin at 15 arc 
minutes. This second-to-largest scale was chosen for this normalization since the largest bin 
has a significant dependence on the second integration constant. 

The B correlator is systematically positive, indicating that some residual non-lensing 
correlations may be contributing. The B mode is about three times smaller than the E 
mode, so this systematic is clearly significantly weaker than the lensing signal. Due to 
symmetry, the E and B error bars are equal in magnitude, with the latter appearing larger 
due to the logarithmic scale. However, the amplitude of the £?-mode is larger than the 
intrinsic alignment prediction which should not exceed a few percents of the lensing signal 
here. Therefore we conclude that most of the 5-mode measured here is due to Point Spread 
Function and/or image analysis residual error. 
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Fig. 1. — Correlation functions. The top panel shows the E-type correlator (boxes) and 
the .B-type correlation (crosses). The bottom panel shows the raw observables: the total 
correlation £ + (boxes) and the differential correlation (crosses). The error bars in the 
bottom panels are fully uncorrelated, while they are mostly uncorrelated in the top panel. 



-7- 




2 
bin separation 



Fig. 2. — Cross correlation coefficient. The two panels show the corresponding cross correla- 
tion coefficient between bins in figure 3. The boxes are the cross correlations for bins of the 
same type (i.e. E — EotB — B), and the crosses between types (e.g. E — B). All covariances 
measured relative to the bin at 4 arcminutes. Other covariances are similar. Not included is 
the integration constant which can move one global offset from all E bins to the B bins. 
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The bin cross correlations are shown in figure 2. The cross correlations are quite small, 
always less than 0.2, except for the auto-correlation of each point, which by definition is 1. 
We have defined the cross correlation coefficient 

°ij = , ^ • (10) 
The E and B correlations are slightly anti-correlated. 



3. Windowed Variance 

From the two correlation functions and we can compute the variances of the 
shear field smoothed with various windows. Van Waerbeke et al. (2001) have measured the 
variances directly from the smoothed maps. Working on the maps is limited by geometry: 
they have many defects which obliges us to mask a fair fraction (~ 20%) of the observed 
area. These masks create many boundaries over the field, which intersect quite often a 
smoothing window, enabling possible edge effects. Integrating over the correlation function, 
this problem does not exist. Crittenden et al. (2000) have presented various strategies to 
derive these variances from correlation functions. 



3.1. Top-Hat Filtering 



The simplest statistic is the top-hat variance: one averages the shear over a disk with 
some radius R, squares it and removes the diagonal terms, such that we effectively measure 
a r.m.s. excess with respect to random alignment of galaxies. We write the smoothed shear 
field as 

7?(*)=/ l t (x-x')d 2 x'. (11) 

J\x'\<R 

The expectation value of the variance of this smoothed field (which is precisely the r.m.s. 
excess of galaxy ellipticities) can be expressed in terms of the shear two point correlation 
function £+(r) 



<(7f) 2 > 



nR A 



2R 



rdr 



2R 2 coa-H—) 
K 2R J 



(12) 



The top-hat variances published by previous authors was the sum of the E and B type. 
Here, again, we can decompose the correlation function £ + = ^ E + and obtain two top- 
hat variance statistics. We show the results in the top panel of Figure 3. The last two 
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Fig. 3. — Windowed Variances. The top panel shows the top-hat variances of the E-mode 
(boxes) and the B- mode (crosses). The error bars are highly correlated, with neighbor- 
ing bins having a 90% cross-correlation coefficient. The bottom panel shows the aperture 
map variances, again with E-mode (boxes) and 5-mode (crosses). The bins are modestly 
correlated, with neighboring cross-correlation coefficient of 50%. 
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B bins are negative, and we took the absolute value to fit on this plot. The error bars are 
computed in the same way as in section 2. The covariance relative to the bin at 6 arcminutes 
is shown in Figure 4. We see that the bins are highly correlated, and there are only about two 
"independent" degrees of freedom when summing over all radii. The B mode is a factor of 3 
smaller than the E mode, again suggesting that the lensing signal is real. This decomposition 
does depend on the integration constant £, which we have chosen to have zero B mode at 
19 arcmin. The interpretation here would be that the data is almost consistent with no B 
mode, but does not rule out a large constant value which cancels a large constant E mode. 
We discuss a statistic below that does not have this degeneracy. 



3.2. Compensated Filtering 

The top-hat variance in (12) is an integral measure of the lensing power, and averages 
over all scales up to the smoothing scale. It has relatively small error bars since it bins 
together the signal on many different scales. It is thus not very sensitive to measure scale 
dependence. While the total top-hat power depends only on galaxy pairs at separations up 
to twice the smoothing radius (i.e. one smoothing diameter), the separation into E and B 
modes is non-local, and depends on large scale correlations. 

A purely local decomposition into E and B variances can be achieved using aperture 
mass estimators (Schneider et al. 1998). Locality here means that the decomposition into 
E- and B-type variances on scale R depends only on galaxy pairs whose separation is less 
than 2R. We define a zero mean mass aperture window 

W(r) = ^(l-r 2 )Q-r 2 ) (13) 
for which J" 1 rdrU(r) = 0. The aperture mass 

7* = ^ / d 2 n t (r)Q(r/R) = -L J d 2 TK{v)U{r / R) (14) 

for the tangential shear 7 t corresponding to 71 in the frame of r. It expresses the aperture 
weighted optical depth k in terms of the observable shear variation. The shear window is 
Q(r) —U{r) — -^2 Jq r'dr'U{r'). These can be expressed as local integrals over the correlation 
functions 

(( 7 |) 2 ) = 7r / rdrW(r)C+(r)+TT rdrW(r)£-(r). (15) 
Jo Jo 

W is the 2-dimensional auto-convolution of U with itself, and VV is defined according to 

Crittenden et al. (2000). The B variance arises by using 72 in (14), or changing the sign in 

the second term of (15), which is zero for weak lensing. 
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Fig. 4. — Cross correlation coefficient. The two panels show the corresponding cross correla- 
tion coefficient between bins in figure 3. The boxes are the cross correlations for bins of the 
same type (i.e. E — EotB — B), and the crosses between types (e.g. E — B). All covariances 
measured relative to the bin at 6 arcminutes. Other covariances are similar. The top panel 
refers to the top-hat, and the bottom panel to the aperture mass. 
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We see the results in the lower panel Figure 3. There is no integration constant, since 
the window has zero mean, so the B mode is an independent estimate. For this reason, 
this plot can be directly compared with Figure 4 in Van Waerbeke et al. (2001): the results 
are consistent but the errors are about twice smaller here, which means that the small 
systematic which was measured in Van Waerbeke et al. (2001) is in fact more significant 
than what we expected. This mean that some among of signal cleaning is required when 
estimating the cosmological parameters with such small error bars (this will be done in a 
forthcoming paper), and also that it is necessary to investigate in more details the possible 
source of systematics in the shape measurement process. 



4. Power Spectrum 

We can measure the power spectrum directly from the correlation function (Kamionkowski 
et al. 1998). The Fourier space power spectra trivially project the power into E modes, which 
are aligned with the wave vector, and B modes which are at 45 degrees. We have 

P E (k) = n 9d9[£ + J (k9) + Z-Ji(k9)], P B (k)=n 9d9[£+ J (k9) - f_ J 4 (k9)} (16) 
Jo Jo 

in terms of the Bessel functions J n . We map k ~ / in spherical harmonic space, and in the 
small angle approximation 1(1 + 1)Ci/2tt = k 2 P(k)/2ir. The inverse mapping is given as 

-I /»oo -I /»oo 

f + (0) = — / kdk[P E (k) + P B (k)}J (k9), f_(0) = — / kdk[P E (k) - P B (k)]Ji(k9) 
2vr Jo ' 2tt J 

(17) 



Numerically, the integral of (16) is dominated by angular scales k ~ 1/9, but the noise 
may receive contributions from all scales, especially scales for which there is little data and 
therefore a lot of noise. To minimize these effects, we used a maximum likelihood power 
spectrum inversion. In analogy to (8) we can express (17) as a linear system 

^~=T ij Pf B ±a i (18) 

We re-scale (18) by a diagonal noise matrix N u = <j,i and parametrize P EB as a sum of band 
powers. We used a sum of 8 basis functions which is proportional to 1/k within each band. 
Tij is then a 20000 by 16 matrix, and we obtain an over-determined set of 20000 equations 
of equal weight for 16 unknowns. The least squares solution gives 

p eb = (rt N t NT yi N -i£+- ( 19 ) 

The correlation function at zero lag is not observable, so the determined angular power 
spectra are indeterminate up to constant. The intrinsic ellipticity white noise has equal E 
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and B modes, so we are free to subtract the same constant white noise component from both. 
We did so by forcing the B mode to be zero at I ~ 5000. This freedom of integration is not 
related to the integration constant in the E-B separation (5), which is one at large scales. In 
the latter case, we chose the integration from infinity rather than integrating from zero. 

The non-observability of the correlation function at zero lag is fundamental. £(0) = 
(e 2 ) + (7 2 ), but we cannot measure the two terms separately. In practice, the first term, 
corresponding to intrinsic ellipticies, is much larger than the rms gravitational shear described 
by the second term. It is perhaps amusing that the weak lensing power spectrum has an 
unobservable offset, while the raw correlation functions do not, in contrast to galaxy 
power spectra where the converse is true. One does obtain lower bounds on the offset since 
power spectra cannot be negative. We chose the bin at I ~ 5000 since the B-mode had 
the smallest relative error bar. This subtraction leaves the higher / slightly negative, but 
consistent with zero within the error bars of each bin. The indeterminacy of the zero lag 
correlation function also manifests itself as an arbitrary integration constant c/r 2 falling off 
as inverse separation squared in the aperture mass statistic, i.e. a line of slope -2 in figure 3. 

The result is shown in Figure 5. For the bin at 4.5arcmin, the _B-mode is consistent 
with zero. At large and small scales, we have a larger amount of systematics. This error 
analysis does not contain cosmic variance, so this plot is not appropriate for measuring the 
slope of the power spectrum. The covariance between bins is shown in Figure 6, showing a 
small covariance. 

We introduce a mixing between modes since the integral (16) is truncated at both large 
and small radii. This coupling is not included in the error covariance. To quantify this 
effect, we generated a pure power law correlation function with the same binning. We used 
£+(#) = oc 1/0, which is a pure E mode (Crittenden et al. 2000). This test generated 

a small negative 5-mode which was largest at the leftmost bin, which was a factor of 17 
smaller in amplitude than the .E-mode, and thus insignificant compared to the other sources 
of error. 

To compare the power spectrum analysis with the aperture mass window from the 
previous section, we have overplotted the square of the fourier transform of the aperture 
window (13) for a radius of 5 arcmin in figure 5 as a dashed line, with the window units 
on the right box side. This window peaks at a scale corresponding to waves slightly longer 
than 5 arc minutes. The window has a minimum at 9 m = ±a/2/36> ~ 4' which has a spacing 
of 8' between the two minima, comparable to a wave of wavelength 8'. We can see the 
qualitative drop of the aperture mass near 20' in Figure 3 as the corresponding drop in the 
power spectrum, where the E mode drops to the level of the B mode. 
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Fig. 5. — Spherical harmonic power spectrum. Again, the boxes indicate the E type power 
spectrum, while the crosses denote the B type power. The angular scale is the spacing of 
peaks of the wave functions at the equator for m=l. The dashed line is the window function 
for the aperture mass statistic at 5', with amplitude labeled on the right bounding box edge. 
The aperture mass is analogous to the integral of the window with the power spectrum using 
equal weight per natural logarithmic interval. 
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Fig. 6. — Bin-bin cross correlation of the angular power spectrum in Figure 5. The covari- 
ances are relative to the bin at 1=1200. 
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The optimal least squares inversion here has ignored the cosmic variance, and is optimal 
in the limit that the signal to noise is very small. For the B mode, this is true at all scales, but 
for the E mode, this procedure is not entirely optimal. A proper power spectrum estimation 
can be done by applying a full maximum likelihood estimate to the data, which will be 
presented in a future paper. 

5. Conclusions 

We have derived several statistics from the two correlation functions measured in the 
VIRMOS-Descart survey. For each statistic, we are able to decompose it into orthogonal E 
and B channels. Since the two channels add to the total power used in previous studies, 
the error in each is ~ \[2 smaller than in the total power. The weak lensing signal is 
solely contained in the E channel, and the B channel provides a monitor to watch for any 
potential systematic contamination. We have analyzed three categories of statistics: the 
correlation function, windowed variances, and angular power spectrum. In each case the 
£>-mode was small, but statistically significant, given the dramatic reduction in the error 
bars compared to the previous analysis. The power spectrum analysis suggested that most 
of the systematics are coming from a constant value C\ which corresponds to a Poisson noise 
contribution. The origin of the systematics is still unclear, but here we developed the tools 
useful for a more detailed analysis, when the amount of data will be large enough. Given that 
the E, B correlation functions and the top-hat smoothed variances are crucially dependent 
on the choice of an arbitrary integration constant, the most robust statistic appears to 
be the aperture mass. Our analysis shows for the first time a direct measurement of the 
projected dark matter power spectrum, which opens a new window, not only for measuring 
the cosmological parameters, but also for a future direct measure of the three-dimensional 
mass power spectrum at small and intermediate scales. 
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Distribution of Dark Matter" of the European Community under contract No. ERBFMRX- 
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